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Abstract 

Rapidity gaps between two hard jets at the Tevatron have been interpreted 
as being due to the exchange of two gluons which are in an overall color- 
singlet state. We show that this simple picture involves unitarity violating 
amplitudes. Unitarizing the gluon exchange amplitude leads to qualitatively 
different predictions for the fraction of t-channel color singlet exchange events 
in forward qq, qg or gg scattering, which better fit Tevatron data. 



I. INTRODUCTION 



Over the past few years rapidity gaps, i.e. pseudorapidity regions without hadronic 
activity, have been observed in hadronic collisions at both the HERA ep collider [1| and in 
pp collisions at the Fermilab Tevatron @-§|]. Such rapidity gaps are widely attributed to 
the exchange of color singlet quanta between incident partons |6|-|9|], the exchange of two 
gluons in a color singlet state being the simplest such model [0 . At the Tevatron, a fraction 
fga P ~ 1% of all dijet events with jet transverse energies E T j > 20 GeV and jet separations 
of more than three units of pseudorapidity exhibit rapidity gaps between the jets. This 
observation is particularly striking since it demonstrates that color singlet exchange effects 
in QCD events are relevant at momentum transfers of order 1,000 GeV 2 , raising the hope 
that perturbative methods can be used for quantitative descriptions. 

A gap fraction of order one percent was in fact predicted by Bjorken [0] , in terms of a 
fraction f s w 0.15 of dijet events which are due to t-channel color-singlet exchange and a 
survival probability P$ of rapidity gaps of order 10% IjTyTOf, 

fga P — f s Ps ■ (1) 

Here the survival probability estimates the fraction of hard dijet events without an underly- 
ing event, i.e. without soft interactions between the other partons in the scattering hadrons. 
Such multiple interactions would fill the rapidity gap produced in the hard scattering pro- 
cess. For Qq elastic scattering, Bjorken estimated the color-singlet fraction f s in terms of the 
imaginary part of the two-gluon t-channel exchange amplitude, which is known to dominate 
the forward scattering amplitude for t-channel color-singlet exchange. In impact parameter 
space, at impact parameters small compared to R = 0(1/ A), the result is 
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Here 2/9 is the relative color factor of the two-gluon color-singlet to the one-gluon color- 
octet exchange cross section and R is an infrared cutoff parameter which regularizes the 



two-gluon loop-integral. This model for the calulation of the color singlet fraction f s , with 
the two gluon-exchange amplitude replaced by its imaginary part, will be called the two- 
gluon exchange model in the following. 

In this model, the color singlet fraction grows with the color charge of the scattered 
partons. For qg and gg elastic scattering f s would be larger by factors 9/4 and (9/4) 2 , 
respectively [f|. This results in a substantial decrease of the observable gap fraction as 
the contribution from gluon induced dijet events is reduced, e.g. by increasing the average 
transverse momentum of the observed jets and thereby the Feynman-x values of the incident 
partons. Such measurements have recently been reported by both the CDF || and the DO 0] 
Collaborations, and no such effect is observed. In fact, the DO data are compatible with a 
slight increase of the gap fraction with increasing jet E T , casting doubt on the validity of 



the two-gluon exchange model |12 



In this paper we reconsider the basic ideas behind the two-gluon exchange model. We 
demonstrate its limitations and show that, even when starting from this perturbative pic- 
ture of rapidity gap formation, the determination of the color singlet exchange fraction f s is 
essentially nonperturbative. We start from a basic feature of the two-gluon exchange model: 
unitarity fixes the imaginary part of the t-channel two-gluon exchange amplitude in terms of 
the Born amplitude and this imaginary part dominates t-channel color singlet exchange . 
Rewriting this relationship in terms of phase shifts, the one- and two-gluon exchange am- 
plitudes are found to be too large to be compatible with unitarity. Phase shift unitarization 
leads to a more realistic description, in which the total differential cross section remains 
unchanged compared to the Born result, but with t-channel color singlet exchange fractions 
which differ substantially from the expectations of the two-gluon exchange model. These 
features are demonstrated analytically for fixed values of the strong coupling constant, a s , 
in Section [TI|. In Section |TTTj we then perform a numerical analysis for running a s , showing 
that the key properties of the fixed-a s results remain unchanged. 

The predicted color singlet fractions are found to very strongly depend on the regulariza- 
tion of gluon exchange at small momentum transfer, however, and thus cannot be reliably 



calculated within perturbation theory. Within our unitarized model the non-perturbative 
effects can be summarized in terms of two parameters, the survival probability of gaps, 
P s , and a universal Coulomb phase shift, ipo. Implications for the formation of gaps at 
the Tevatron are analyzed in Section [TV[ In particular we calculate how the gap fraction 
between two hard jets varies with jet transverse energies and jet pseudorapidity separation 
and then compare predicted fractions with Tevatron data f4j||. Our conclusions are given 
in Section |V|. 

II. ELASTIC SCATTERING AMPLITUDE AND UNITARIZATION 

Consider the elastic scattering of two arbitrary partons, p and P, 

p{ii) + PUi)^P{i2)+P{32) , (3) 

at momentum transfer Q 2 = —t. Here . . . , j'2 denote the colors of the initial and final 
state partons. The cross section and the partial wave amplitudes are completely dominated 
by the forward region, Q 2 <C s, where the Rutherford scattering amplitude, 

M = -87ta s S -T a ® T' a = 87ra s ^ F c = M F c , (4) 

provides an excellent approximation. Note that helicity is conserved in forward scattering, 
hence spin need not be considered in the following. The only process dependence arises from 
the color factor F c = T a ® T' a . 



A. Diagonalization in Impact Parameter and Color Space 

In order to study unitarity constraints, we need to diagonalize the amplitude in both 
momentum/coordinate space and in color space. The first step is most easily achieved by 
transforming to impact parameter space, 



Neglecting multi-parton production processes, i.e. inelastic channels, unitarity of the S- 
matrix implies the relation 

JmT(b) = ^|T(b)| 2 , (6) 

for the full 2^2 scattering amplitude T(b). Eq. (||) represents a matrix relation in color 
space. More fully it can be written as 

2m T(h) i 2 j 2 ,iiji — 4^ 

Tt(b) (7) 

where the sum runs over the dimension of the color space, dc — 9 for Qq and Qg scattering 
and d c = 24 (64) for qg (gg) elastic scattering. 

Since the color factors can be written as hermitian matrices, the right-hand side of Eq. (|7|) 
represents a simple matrix product of the color matrices. This product is easily diagonalized 
by decomposing the color factors F c into a linear combination of projection operators onto 
the irreducible color representations which are accessible in the s-channel, 

F c = {Fc)i 2 j2,hji = £ fk (Pk)i 2 j2,iiji = fk Pk ■ (8) 
k k 

For the case of quark- ant iquark elastic scattering, for example, with color decomposition 
3 (g) 3 = 1 © 8, the color factor can be written in terms of Gell-Mann matrices as 



'A a \ /A a \ 4 1 /A a \ /A a " i i 

{Fc)i 2 j2,iiji = ( ~2~ ) i2il ( ~2~ J 5 ' 1J ' 2 = g^ii^i^wa ~ 3 \ ~2~ / \ ~2~ / i2 ^ 2 = 3 ^ l ~ 6 ^ 8 



For all cases, Qg, Qg, gg and elastic scattering, the decomposition into s-channel projec- 
tors is summarized in Table |. This color decomposition, combined with the transformation 
to impact parameter space, diagonalizes the unitarity relation for elastic scattering ampli- 
tudes. 



B. Phase Shift Analysis at Leading Order in a s 



Expanding the full 2^2 amplitude T(b) into s-channel projectors, T(b) = J2k ^fc(b) Pk, 
the individual coefficients are seen to satisfy the unitarity relation (|6|). The full scattering 
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TABLE I. Representations and color operators for QCD elastic scattering. The indices of the 
projection operators Pk in the last column represent the dimensionalities of the irreducible color 



representations in the s-channel. Results for the 8 (g> 8 decomposition are taken from Ref. [13]. 
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amplitude can thus be written in terms of real phase shifts Sk{h) for specific color represen- 
tations in the s-channel, 



T(b)=^T fc (b)P fc = -^2^ (e 



2«<5 fc (b) 
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(10) 
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Within perturbation theory, the individual phase shifts Sk{h) can be expanded in a power 
series in a s . The lowest order term is fixed by the Fourier transform of the Born amplitude 
(|), which, however, diverges at small |q| and needs to be regularized. This is most easily 
done by an infrared cutoff, |q| > 1/-R, of the integral ([Sp. One can interpret this infrared 
cutoff as a consequence of confinement; the color singlet nature of hadrons at scales larger 
than pd 1 fm does not allow long wave-length gluons to couple and, hence, soft gluon exchange 



must be suppressed. The cutoff R is related to the size of the hadronic wave-function |TJ 
and can be considered as a nonperturbative parameter in the following. 

With the cutoff |q| > 1/R, the Fourier transform of the Rutherford amplitude M.® to 



impact parameter space is given by 

T (b) = 4 S ^(bg^ + 2(log2- 7 ))F c = As6 (b)F c . (11) 

Here 7 = 0.577215 ... is Euler's constant, and terms of order h 2 /R 2 are neglected. Compar- 
ison with ( pup yields 

5 k (b) = f k 5 (h) + 0(a 2 ), (12) 



where the f k are taken from Table |. Keeping the lowest order term in fll2|) only, the 
transformation back to momentum space can be performed analytically for the full amplitude 
in flT0|), with the result 

M(q) = 8n as ^Y,hPkexp(t as f k logR 2 q 2 + 0((a s f k ) 3 )) . (13) 

1 k 



As in the analogous QED case JTH], the coefficient M. k of each projector P k is just the 



Born amplitude, multiplied by an infrared divergent phase factor i.e. 

M=J2 M k P k = fkM e* M P k , (14) 

k k 

with 

ii = a s \ogR 2 q 2 + 0(a 3 s ) . (15) 

This general structure has important consequences for the total differential cross section, 
summed over all colors, and for the t-channel color singlet exchange rate. From ( ITS ) the 
color summed amplitude squared is given by 

£ \M\ 2 = Y,\M k \ 2 d k (16) 

colors k 

where d k is the dimensionality of the fcth irreducible color representation in the s-channel 
(see Table |). Since each M. k equals its tree level value f k A4 , up to a phase, the total 
differential cross section remains unchanged by our phase shift unitarizationQ 



*This is not necessarily true for other unitarization prescriptions. For example, for a "iT-matrix 
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C. Color Singlet Exchange Fraction 

In order to understand the rapidity gap rate in hadronic collisions, we need the fraction 
of dijet events which are produced without color exchange in the t-channel. The color factor 
describing this situation is given by 

F s = 4m <W = 1 ? ( 17 ) 

which is just the unit operator, as far as the decomposition into s-channel projectors is 
concerned. For any given 2 — > 2 process we define the t-channel color singlet exchange am- 
plitude, Ai s , as the coefficient of F s . The exchange of color octet quanta in the t-channel or 
of yet higher color representations will be orthogonal to this term, i. e. no interference terms 
arise once the squared amplitude is summed over all colors, and this makes the definition of 
M. s unique. Decomposing the full amplitude as M. = J2k-Mk Pk = -Ms F s + . . ., the color 
singlet fraction, f s of Eq. (0), for any particular 2 — > 2 process, is then given by 

n _ Scolors | J2kMk d k \ 2 , s 

h ~ Ecolorsl-MI 2 ~ d C Ek\M k \ 2 d k ' 1 } 

where dc = Sfc d k is the dimensionality of the full color space. 

Let us apply this expression to the full, unitarized amplitude of Eq. (|T4|). The lowest order 
Rutherford amplitude M.$ cancels in the ratio of Eq. (|TE|), which hence can be evaluated in 
terms of the coefficients fk and the dimensionalities dk of Table |. We find 

UQq^Qq) = j| sin2 ^ (19) 
fs(Qq^Qq) = f sin 2 ^, (20) 



unitarization" , with ?fc(b) = 4s/fcTo(b)/(4s — ifhT^ih)), we numerically find substantially reduced 
unitarized cross sections. For gg scattering, da(gg — * gg) / dcosO can be a factor 5 lower than the 
Born result, even for large scattering angles. Such a unitarization procedure would be completely 
unacceptable phenomenologically. 
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fs(qg ^qg) = - sin 2 - + — sin 2 ip , (21) 

9 2 1 2 3^ 82, 

tAqq — > go) = — sin — sin — H sin 2w . (22) 

jsyyy yyj 16 4 16 4 128 ^ V ' 

These color singlet fractions are plotted as a function of ip in Fig. [I]. 




FIG. 1. Fraction of i-channel color singlet exchange events in Qq, Qq, qg and gg elastic 
scattering as a function of the universal phase tp. See text for details. 

From (|15|) we see that, formally, the universal phase ip is of order a s and thus the color 
singlet fractions are of (9(a 2 ), which agrees with Bjorken's result of Eq. (|2|). Expanding ( |20|) 
to lowest order, one obtains 

fs(Qq^Qq) « \ i? = I |cUogi? 2 q 2 | 2 , (23) 

which appears to be four times larger than the result given in (fj). This conundrum is resolved 
by observing that (Q) represents a cross section ratio in impact parameter space while (|T9|- 
|2"2"D are the color singlet fractions in momentum space. Indeed, in impact parameter space 
and to leading order in a s , Eqs. (|ll3| - |i"2"D imply 

9 



f^{Qq^Qq) « 2 - |5 (b)| 2 = 2 - 



a s . R 2 |2 

T log b^ 



(24) 



which agrees with Eq. (0) .[] The factor four difference between the color singlet fractions in 
momentum and impact parameter space can be traced to a binomial factor 2 in the Fourier 
transform to momentum space of log 2 R 2 /b 2 = \og 2 (qR) 2 — 21og(6g) 2 log(g-R) 2 + log 2 (6g) 2 
compared to the transform of log R 2 /b 2 = log(qR) 2 — log(6g) 2 : the relevant term is the one 
linear in log(6g) 2 , which is enhanced by a factor 2 in the first case. Since all experiments 
are analyzed in momentum space, the predicted color singlet fractions in the two-gluon 
exchange model need to be multiplied by a factor 4 compared to the results derived from 
(0). For gluon-gluon scattering in particular, this would lead to a color singlet fraction of 
f s (gg) = 0.15 x 4 x (9/4) 2 « 3, which obviously cannot hold. 

The reason for this problem is the fact that, with the same arguments as used in Eq. (0), 
the phase appearing in Eqs. (|19|-|2"2"|) is approximately given by i\) ~ 127r/(33 — 2rif) = 1.64 



which is too large to make use of a small angle expansion. Instead, ( |22"D predicts f s (gg — > 
gg) = 0.39, which is still a surprisingly large color singlet fraction but almost an order of 
magnitude smaller than the two-gluon exchange result. 

Indeed, a resummation of higher order terms is required by unitarity. In the two-gluon 
exchange model, in impact parameter space, the amplitude for a particular color represen- 
tation in the s-channel is given by 

^r fc (b) = S (b)f k + i(5 (b)f k ) 2 (25) 
4s 

with So(h) ~ 67r/(33 — 2nf) = 0.82, for small impact parameters. Here box corrections to 
the real part are neglected. The unitarity relation (H), on the other hand, implies 



7& ( j-Ub) 



< 0.5 , (26) 



^For Qq — ► Qq scattering and in impact parameter space we thus agree with Ref. while the 
additional factor four in momentum space was missed in Ref. 0. 
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a condition which is violated for all color channels with \fj.\ > 2/3 in Table |. Since the two- 
gluon exchange model violates unitarity, we study implications of the unitarized extension 
of this model in the following. 

So far we have estimated the value of the universal phase ip by replacing a s by its running 
value a s (Q 2 ) in (JTSp, 

m) - «.«« H w = ^ (i + = | + a . ( ^) io g «v , ( 2 7) 

and then using the asymptotic expression for logQ 2 — > oo, ie. neglecting the 0(a s ) term 
in fl2"7|). This corresponds to setting the cutoff parameter R = 1/A. In a more complete 
calculation i? describes the transverse length scale at which color screening, due to other 
partons in the proton, sets in, thus suppressing the effective gluon coupling. In effect, 
k = RA « 1 is a non-perturbative parameter for which we only have a rough guess, and 
which may be uncertain to at least a factor three, if not an order of magnitude. With 
4"7r/ ' (5q = 0.527T = 1.64, and a s (Q 2 ) ~ 0.14 at momentum transfers relevant for the Tevatron 
rapidity gap data, a variation of k by a factor 10 leads to changes in if) by 30% or more. 
A change of this order, in particular an increase of if>, can drastically change the predicted 
color singlet fractions for individual processes, as is obvious from Fig. |l|. 

III. RUNNING COUPLING EFFECTS 

An increase of if) is, in fact, to be expected when including running coupling effects in the 
determination of the tree level phase shift 5o(b) in Eq. (|TT|). A running coupling increases 
the average size of the Born amplitude ([|) in the Fourier transform to impact parameter 
space, which leads to larger values of <5o(b) in (0) and this translates into a larger phase if) 
in©. 

We have analyzed this question quantitatively by determining the partial wave phase 
shifts which correspond to the Born amplitude, with running a s (Q 2 ), and then unitarizing 
the partial wave amplitudes as in flTUp. In terms of the Born amplitude A4q(Q 2 ) the phase 
shifts for fixed angular momentum J are, to lowest order, given by 
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6j=—Jdcos9Pj(cos9)M (Q 2 = ~(l-cos9)) . (28) 

This integral is singular at Q 2 = 0, via the 1/Q 2 pole of the Rutherford amplitude, and 
at Q 2 = A 2 , via the Landau pole of a s (Q 2 ). Both singularities need to be regularized, for 
which we introduce two (independent) mass parameters, M and M a . We thus replace the 
Born amplitude by 

AVQ 2 ) = 8. a ,(Q\ AfJ ^ = j-^ ^ . (29) 

We do not expect this amplitude to correctly describe the actual QCD matrix elements 
at low Q 2 . By varying the infrared cutoff parameters M and M a we rather explore the 
importance of the small Q 2 region and, thus, the importance of non-perturbative effects 
which we are unable to calculate. The resulting unitarized amplitudes are now given by 

oo 

M k (6) = 8tt ]T (2 J + l)Pj(cos 6)%,j (30) 
j=o 

with 

/ 4t f 1 1 \ 

%T KJ = exp li —f k / dxPj(x)- — - 1 . (31) 

V Po J -i (z - x) log {z a -x)J 

Here the Pj(x) are Legendre polynomials, and z = 1 + 2M 2 /s and z a = 1 + 2M 2 /s contain 

the two regularization parameters. 

The integration in (|3l|) and the partial wave sum in (|30D are performed numerically. 
Results are shown in Fig. [|, where possible variations of the color singlet fraction with the 
regularization parameters M and M a are explored. A very strong variation of f s is found for 
individual scattering processes. But, as we shall demonstrate below, these variations tend 
to be averaged out to a large extent when summing over the various partonic subprocesses 
which contribute to actual dijet data. 

The parameter dependence of f s for running strong coupling constant in Fig. ^ is rem- 
iniscent of the variation of f s with the universal phase ip, which is shown in Fig. [I] for our 
analytical results (|TU|-P^). The similarities between the running coupling partial wave ex- 
pansion (RPWE) and the impact parameter space (IP) results go far deeper, in fact |16 
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FIG. 2. Dependence of the i-channel color singlet-exchange fraction on the regularization 
parameters (a) M and (b) M a (see Eq. (^)). Results are shown for fixed momentum transfer, 
Q = 30 GeV, and fixed forward scattering angle. 

Similar to the analytical IP results, the RPWE calculation leads to a total differential cross 
section, summed over all colors, which agrees with the Born result to few percent accuracy. 
In addition, the color singlet fractions at fixed Q 2 are found to be independent of scattering 
angle or parton center of mass energy. Indeed, with an accuracy of a few percent, the nu- 
merical results of the RPWE calculation can be parameterized as in (HI-ID, with a phase 
■0 = ip(Q 2 , M, M a ) which is universal for all subprocesses. As is obvious from Fig. |2| the 
dependence of ip on the regularization parameters is quite strong. Its Q 2 -dependence, on 
the other hand, is logarithmic only and the ansatz 



0(Q 2 ) = tfj + ih log % + 2 log 2 % 



(32) 



Qo Qo 

provides an excellent parameterization of the RPWE results. Representative values for 
the coefficients ipi are given in Table [H]. Only ijj is found to depend appreciably on the 
regularization parameters while ip\ = 0.256 and ip2 — —0.019 (for Q = 50 GeV) are 
constant within the numerical uncertainty. The variation of ip(Q) with momentum transfer 
is modest, with Aip(Q) « ±0.2 in the interval 20 GeV < Q < 100 GeV. In this region, which 
is of interest for the Tevatron, the RPWE values for ip(Q 2 ) in general are substantially larger 
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TABLE II. Dependence of the universal phase ip(Q 2 ) on the regularization parameters M and 



M a (see Eq. (p9|)). The ipi are the coefficients of the expansion in (|32| ) for Qq = 50 GeV. 
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than for the analytic impact parameter space calculation, thus confirming the qualitative 
arguments made earlier. 

Beyond demonstrating the relation of the numerical RPWE results to the analytical 
expressions derived in impact parameter space, the existence of the simple parameterization 
in fl32|) is very important in order to compare our calculations to experimental data. The 



numerical integrations which need to be done in (|31|) are too slow to be performed for 
individual phase space points in a Monte Carlo calculation of dijet cross sections. With the 
above observations this is not necessary, however, since instead we can use the analytical 
results of ( |T§| - P2D together with the parameterization of Eq. 
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IV. COMPARISON WITH TEVATRON DATA 



Both the DO Q and the CDF [H] collaborations at the Tevatron have analyzed the fraction 
of dijet events with rapidity gaps, as a function of both the transverse energy, E T , and the 
pseudorapidity separation, Arj = \rjj 1 — rjj 2 \, of the two jets. As these phase space variables 
change, the composition of dijet events varies, from mostly gluon initiated processes at 
small Et and At] (and, hence, small Feynman-x, xf) to Qq scattering at large values. A 
dependence of the gap fraction on the color structure of the scattering partons would thus 
be reflected in a variation with Et and At]. 

Bjorken's two gluon exchange model, which is equivalent to the small ip region in our 
analysis, predicts a larger fraction of color singlet exchange events for gluon initiated pro- 
cesses (see Fig. [I] for ip < 0.37r). The gap fraction should thus decrease with increasing 
Et or Arj. The opposite behavior is expected in statistical models of color rearrange- 
ment pT| , [12| . Here the eight color degrees of freedom for gluons, as compared to three for 
quarks, make it less likely for gluon initiated processes that t-channel color singlet exchange 
is achieved by random color rearrangement. This would lead to a smaller gap fraction at 
small xp and therefore small Et or Ar]. 

In the unitarized RPWE framework, the dependence on the regularization parameters is 
sufficiently strong to encompass both scenarios. This is demonstrated in Figs. |^ and ^, where 
the results of the running coupling analysis for three choices of the regularization parameters 
are compared with Tevatron data, taken at ^Js = 1800 GeV. The data correspond to dijet 
events with two opposite hemisphere jets of Et > 20 GeV, \rjj\ > 1.8 (CDF) or Et > 30 GeV, 
\rjj\ > 1.7 (DO). DO data are taken from Ref. 0] and show the fraction of dijet events with 
rapidity gaps. CDF || shows the ratio of gap fractions in individual Arj and Et bins to 
the overall gap fraction in the acceptance region. For comparing our calculation with the 
data we fix the survival probability P s in Eq. (P to reproduce the overall gap fraction in the 
acceptance region, which was measured as f gap = (0.85±0.06±0.07) % for the DO sample and 
fgap — (1-13 ± 0.12 ± 0.11)% for the CDF sample. Required survival probabilities strongly 
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depend on regularization parameters and vary between 1.9% and 5.5%, which is on the low 
side of previous estimates [f^JTOfl . For a given choice of regularization parameters, predictions 
for the Et or At] dependence of the gap fraction are quite similar for the CDF and DO cuts. 
To the extent that the two data sets are consistent within errors, it is not yet possible to 
discriminate between different choices of regularization parameters, i.e. to obtain sensitivity 
to the non-perturbative dynamics. 

DO data somewhat favor color singlet fractions which grow with xf and which are more 
in line with expectations from color evaporation models []lT , |l2 |. Note that our unitarized 
gluon exchange model, with M = 0.2 GeV and M a = 0.5 GeV is able to describe this trend, 
even though it is an extension of the two-gluon exchange model. CDF data slightly prefer a 
gap fraction which decreases with increasing At] and, hence, with larger xp. The unitarized 
gluon exchange model, with M = 0.1 GeV and M a = 0.2 GeV describes such a situation. 
Comparison with Fig. shows that this set of parameters predicts a much smaller gap 
fraction for Qq scattering than for gluon initiated processes, which is qualitatively similar 
to the two-gluon exchange model Indeed, the shape of the gap fraction for the two 

gluon exchange model is very similar to the long-dashed curves in Fig. [|. 

Clearly, the data are not yet precise enough to unambiguously distinguish between these 
different scenarios. On the theoretical side, the variation of the RPWE predictions with 
model parameters highlights the limitations of a perturbative approach to the color singlet 
exchange probability. Taking the phase ipo and the survival probability P s as free parameters, 
the unitarized two-gluon exchange model is clearly capable of fitting the present Tevatron 
data, however. 



V. DISCUSSION AND CONCLUSIONS 

The formation of rapidity gaps in hadronic scattering events is a common occurrence, and 
its ubiquity asks for a theoretical explanation within QCD. The formation of gaps between 
two hard jets at the Tevatron is particularly intriguing and is commonly being explained in 
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FIG. 3. Dependence of rapidity gap fraction as a function of (a) Arj and (b) jet Et for three 
choices of regularization parameters which are given in units of GeV. Symbols with error bars 
correspond to the DO data as given in Ref . [||] . Also shown are the survival probabilities needed to 
reproduce the overall rapidity gap fraction of 0.85% as measured by DO. 

terms of color singlet exchange in the t-channel, be it via an effective color singlet object 
like the "Pomeron" or via a statistical color rearrangement, in terms of multiple soft gluon 
exchange. 

"Pomeron" exchange models build on the observation that color singlet exchange in 
the t-channel can be achieved in QCD via the exchange of two gluons, with compensating 
colors || . When trying to build a quantitative model for the formation of rapidity gaps [|7| , 
one encounters infrared divergences in the color singlet hard scattering amplitude, which 
in a full treatment would be regularized by the finite size and the color singlet nature of 



physical hadrons ||14j| . In turn, this indicates that non-perturbative information may be 
indispensable for a quantitative understanding of the hard color singlet exchange process. 

We have analyzed this question within a particular model, based on the unitarization 
of single gluon exchange in the t-channel. The Low-Nussinov model || corresponds to a 
truncation of the unitarization at order a 2 s . We find that, for any reasonable range of regu- 
larization parameters, the two-gluon exchange approximation violates partial wave unitarity, 



17 



o 

CO 

~o 

CD 
N 



O 

c 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 




Vs=1800 GeV E T >20GeV 1.8<|tiI<3.5 

p' =2.88'% M=0.1 M '=0.2 ' 
P s =5.51 % M=0.2 M a =0.5 
P s =2.46% M=1.0M a =2.0 



Vs=1800 GeV E T >20GeV 1.8<|tiI<3.5 



----- P s =2.88 % k=0A M =6.2 
P s =5.51 % M=0.2M a =0.5 
P s =2.46% M=1.0M a =2.0 



4.5 



5 5.5 

ATi=|Ti r Ti 2 | 



6.5 



20 25 30 35 40 45 50 55 60 
E T [GeV] 



FIG. 4. Dependence of rapidity gap fraction as a function of a) Ar] and b) Et- Shown is 
the gap fraction in a given bin, normalized to the inclusive gap fraction of 1.13% as measured by 
CDF. CDF data Q are shown with error bars together with predictions of the unitarized RPWE 
calculation, for the three choices of regularization parameters as in Fig. [| 

and thus a fully unitarized amplitude is needed for phenomenological applications. 

The unitarization of hard elastic quark and gluon scattering is not unique, of course, 
but any acceptable method must preserve the successful description of hard dijet events 
by perturbative QCD. The phase shift approach used here fulfills this requirement: the 
unitarization does not change the Born-level predictions for the color averaged differential 
cross sections. As a corollary, the color-inclusive dijet cross section is independent of the 
regularization parameters which need to be introduced for the full phase shift analysis. 

The situation is entirely different when considering the t-channel color singlet exchange 
component which is introduced by the exchange of two or more gluons or by unitarization. 
The t-channel color singlet exchange fraction, f s , is strongly affected by the full unitariza- 
tion and deviates from the expectations of the two-gluon exchange approximation, changing 
even the qualitative predictions of the Low-Nussinov model. These strong unitarization 
effects are reflected by a strong dependence on the precise regularization procedure. This 
cutoff dependence, which parameterizes non-perturbative effects, does not allow to make 
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quantitative predictions for the color singlet exchange fractions in particular partonic sub- 
processes. These limitations, which have been demonstrated here for the two-gluon exchange 
approximation to the Pomeron, may be generic to Pomeron exchange models, and should 
be analyzed more generally. 

In spite of these limitations, we find some intriguing features of the unitarized gluon 
exchange amplitudes. For all partonic subprocesses and for all regularization parameters, 
the color singlet exchange fractions can be described in terms of a single universal phase, 
vp(Q 2 ), which absorbs all non-perturbative effects. This suggests a unified phenomenological 
description of the rapidity gap data, via the gap survival probability P s and the phase ip(Q 2 ). 
Such an analysis goes beyond the transverse momentum and pseudorapidity dependence of 
rapidity gap fractions which have just become available, and should best be performed 
directly by the experimental collaborations. 
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